Shifts in fungal community diversity and potential function under natural forest succession and planted forest restoration in the Kunyu Mountains, East China

Abstract Soil fungi participate in various ecosystem processes and are important factors driving the restoration of degraded forests. However, little is known about the changes in fungal diversity and potential functions under the development of different vegetation types during natural (secondary forest succession) and anthropogenic (reforestation) forest restoration. In this study, we selected typical forest succession sequences (including Pinus densiflora Siebold & Zucc., pine‐broadleaf mixed forest of P. densiflora and Quercus acutissima Carruth., and Q. acutissima), as well as natural secondary deciduous broadleaved mixed forests and planted forests of Robinia pseudoacacia on Kunyu Mountain for analysis. We used ITS rRNA gene sequencing to characterize fungal communities and used the FUNGuild database to predict fungal functional groups. The results showed that forest succession affected fungal β‐diversity, but not the α‐diversity. There was a significant increase in Basidiomycota and a decrease in Ascomycota in the later successional stage, accompanied by an increase in the functional groups of ectomycorrhizal fungi (ECM). Conversely, planted forests exhibited decreased fungal α‐diversity and altered community compositions, characterized by fewer Basidiomycota and more Ascomycota and Mucoromycota. Planted forests led to a decrease in the relative abundances of ECM and an increase in animal pathogens. The TK content was the major factor explaining the distinction in fungal communities among the three successional stages, whereas pH, AP, and NH4 + were the major factors explaining community variations between natural and planted forests. Changes in vegetation types significantly affected the diversity and functional groups of soil fungal communities during forest succession and reforestation, providing key insights for forest ecosystem management in temperate forests.


| INTRODUC TI ON
Soil fungi play a vital role in soil nutrient cycling, are closely related to plants, and regulate plant community dynamics and productivity in forest ecosystems (Harris, 2009;Tomao et al., 2019).They are extremely diverse and functionally versatile and mainly comprise symbiotic, saprotrophic, and pathogenic fungi.Mycorrhizal fungi, such as arbuscular mycorrhizal fungi (AMF) and ectomycorrhizal fungi (ECM), are symbiotic with plant roots and provide nutrients and water for their hosts in exchange for plant carbon (Brundrett & Tedersoo, 2018).Mutualistic symbiotic interactions promote host plant growth (Franco et al., 2014;Roger et al., 2013), improve plant resistance and adaptability to stress (Begum et al., 2019), and create an underground network among plants to transport nutrients (Taylor, 2006).Saprotrophic fungi can degrade soil organic matter (SOM) and plant litter and supply soil nutrients (i.e., nitrogen and phosphorus) to support plant growth in most forest ecosystems (Talbot et al., 2013).Pathogenic fungi generally have negative interactions with plants, inhibiting plant growth and altering plant diversity and community composition (Bagchi et al., 2014).Therefore, soil fungi are of considerable importance in promoting plant growth and supporting ecological stability in forest ecosystems.
Fungal communities are susceptible to various environmental changes in climate (Martíınez-García et al., 2017), land use (Spurgeon et al., 2013), plant species (Dassen et al., 2017;Hedĕnec et al., 2023), and soil physicochemistry (Glassman et al., 2017;Wang et al., 2019).The changes in plant species and soil environment could lead to substantial shifts in the soil fungal community.Changes in plant species during forest succession or reforestation can alter the soil microenvironment by introducing litter and root exudates (Han et al., 2021;Smith et al., 2015).During succession from coniferous forest to broad-leaved forest, broadleaf litter has higher nitrogen and phosphorus content and a lower carbon and nitrogen ratio than coniferous litter.This can lead to a higher litter decomposition rate to facilitate the accumulation of soil-available nutrients (Augusto et al., 2015).Therefore, the available soil nutrient content of broad-leaved forests is sufficiently high to support fungal community growth.Changing plant species composition can result in the recruitment of species-specific microbial groups via the production of root exudates (Santoyo, 2022).Therefore, changes in plant species composition with natural forest succession or reforestation are closely associated with changes in the diversity and community structure of soil microorganisms.Substantial changes in soil fungal species diversity have been observed during plant succession (Ren et al., 2019;Reyes et al., 2019).
However, the varied patterns can be inconsistent, which may be related to the vegetation type and soil environment (Balami et al., 2021;Zhao et al., 2023).Fungal community composition is highly sensitive to changes in successional stages (Bachelot et al., 2018).For example, a universal shift in fungal communities occurs from Ascomycota to Basidiomycota-dominant communities during succession (Jiang et al., 2021;Yan et al., 2020).Changes in the fungal community directly affect its function, which, in turn, drives forest community succession via soil nutrient availability and plant-soil feedbacks.
A high relative abundance of soil ECM in mid-to late-successional stages leads to positive biotic plant-soil feedbacks and strongly promotes the performance of slow-growing coniferous seedlings to favor successional development (Zhao et al., 2023).However, to date, most prior studies have focused on the changes in fungal diversity and community composition during vegetation succession.In this context, relatively little is known about the changes in the functional characteristics of fungi along the forest successional stages.
Understanding the changes in community dynamics and potential functions of fungi during forest succession can provide more direct insights into how changes in the environment affect microbial processes related to forest ecosystem functions.
Reforestation is an important way of forest vegetation restoration.The ecological process of reforestation leads to rapid changes in the aboveground vegetation in terms of tree species, plant diversity, community composition, litter, and soil properties, which play an important role in controlling soil fungal communities.However, the effects of planted forests on soil fungal communities can vary (Carson et al., 2010;Zhang, Li, et al., 2023).For instance, there can be higher fungal diversity and production in planted Pinus forests than in natural Quercus forests (Oria- de-Rueda et al., 2010).Planted forests of Eucalyptus globulus or Pinus pinaster significantly alter the fungal community structure (Carson et al., 2010).Therefore, to date, our understanding of the effects of reforestation (especially planting trees) on soil fungal communities has remained relatively limited.In this context, it is necessary to investigate the responses of soil fun- soil fungal community.Therefore, we hypothesized that soil fungal community diversity and function varied during secondary forest succession and reforestation, and their variation pattern differed.In conclusion, the primary aims of this study were to (1) clarify how the diversity and composition of soil fungal communities change during forest succession and reforestation; (2) explore the changes in fungal functional groups during forest succession and reforestation; and (3) identify the key environmental factors regulating their diversity and functional shifts.

| Study site
The study site is located in the Kunyu Mountain National Nature Reserve (37.20° to 37.31° N, 121.62° to 121.85° E), Shandong Province, China.This area is characterized by a warm temperate monsoon climate, with an average annual temperature of 11.9°C, an average annual precipitation of 984.4 mm, and a frost-free period of 200-220 d.The soil is classified as brown soil.Since the 1950s, this area had suffered from large-scale damage of pine caterpillars, leading to many diseased trees were felled.Then reforestation was carried out by various methods, such as tree and shrub mixed forest, coniferous and broad-leaved mixed forest, and 80% of the original pure pine forest was turned into mixed forest.In the middle and late 1980s, this area was closed to forest cultivation, which accelerated secondary succession forming a majority of P. densiflora forest, interspersed with mixed forests larch and deciduous trees.In this study, five main forest types and non-forest (NF) land were selected.The main forest types include natural pine forests of P. densiflora (PD), natural pine-broadleaf mixed forests of P. densiflora and Q. acutissima (PQMF), and natural broadleaved forests of Q. acutissima (QA), which constitute a natural secondary succession series, and natural secondary deciduous broadleaved mixed forests (MF), including Sorbus alnifolia (Sieb.et Zucc.)K. Koch, Ailanthus altissima (Mill.)Swingle, Kalopanax septemlobus (Thunb.)Koidz., etc., and planted forest of R. pseudoacacia (RP).

| Soil sampling
In November 2020, five field plots (10 m × 10 m) were established within each forest type, and each plot was approximately 50 m apart.
We selected five sampling points in each field plot, taking the center of the square as the sampling center, and four samples equidistant from the central sampling point were collected.Soil samples were collected from the 0 to 10 cm depth of surface soil in each field plot using a stainless-steel collar and were mixed into one composite sample.In total, five soil samples were obtained in five plots for each forest type, thus we collected a total of 30 soil samples (5 forest types × 5 samples +1 non-forest × 5 samples).Each composite soil sample was divided into two subsamples.One sample was immediately frozen in liquid nitrogen, transported to the laboratory as soon as possible, and stored at −80°C until DNA extraction.The other sample was placed in sterile bags and stored at −20°C for the determination of soil physicochemical properties.

| Soil property determination
Soil pH was determined in the supernatant of 1:2.5 soil-water mixtures with a pH meter (Mettler, Switzerland).SOM content was measured following the K 2 Cr 2 O 7 -H 2 SO 4 oxidationreduction colorimetric method.Soil total nitrogen (TN) was determined by a Vario Micro Cube elemental analyzer (Elementar).Total phosphorus (TP) was measured by the molybdenum blue method after wet digestion with H 2 SO 4 + HClO 4 .Total potassium (TK) was determined following digestion with HNO 3 and HClO 4 by a flame photometric detector.Ammonium (NH 4

+
) and nitrate (NO 3 − ) were extracted with 2 M KCl solution using a 1:10 soil/extractant ratio and measured using a nutrient autoanalyzer (Seal, Germany).Available phosphorus (AP) was determined via sodium bicarbonate extraction using the molybdenum blue method.Available potassium (AK) was extracted with ammonium acetate and quantified using a flame photometric detector.All the soil properties were determined according to the method described by Bao (2000).

| DNA extraction, PCR amplification, and high-throughput sequencing
Total genomic DNA was extracted from 0.5 g soil using a FastDNA Spin Kit (MP Biomedical, USA), following the manufacturer's instructions.
A NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific, USA) was used to determine the concentration of the extracted DNA.The ITS1 region of the fungal rRNA gene was amplified using the forward primer ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and the reverse primer ITS2R (5′-GCTGCGTTCTTCATCGATGC-3′).

| Sequence data processing
Raw sequences were trimmed and quality-filtered using Q IIME (v.1.9.1) (Caporaso et al., 2010).Joined paired-end reads were filtered using a series of quality standards, including a quality score higher than 20, length longer than 200 bp, no ambiguous bases, and homopolymers less than 6 bp.The remaining sequences were assigned to operational taxonomic units (OTUs) using Uparse v.11 with 97% sequence similarity (Edgar, 2013), and chimeric sequences were identified and removed.Singleton sequences were excluded from subsequent analyses.The taxonomy was assigned using the UNITE (version 8.0) fungal ITS database (Nilsson et al., 2019).We used FUNGuild (Nguyen et al., 2016) to analyze the functional groups of the fungi in the soil.
Only the two levels of "Highly Probable" and "Probable" were retained for subsequent analyses.Raw data were deposited in the NCBI SRA database under the accession number PRJNA1032641.

| Statistical analyses
One-way analysis of variance and Tukey's HSD test were performed to determine the differences in soil properties, relative abundances of fungal taxa, and functional groups using SPSS (v20.0,SPSS Inc., USA).The α-diversities of the fungi were calculated based on random resampling of 39,600 reads per sample, including the Shannon, Simpson, ACE, and Chao indices.Nonmetric multidimensional scaling (NMDS) was performed to evaluate the variations in soil fungal β-diversity based on the Bray-Curtis dissimilarity matrices at different successional stages and between natural and planted forests.
Analysis of similarity (ANOSIM) was conducted to examine significant differences in the fungal community composition (Clarke, 1993).A distance-based redundancy analysis (db-RDA) was performed to examine the influence of environmental variables on fungal communities and functional groups.Spearman's correlation analysis was performed to explore the associations between α-diversity, or relative abundance of fungi, and soil physicochemical characteristics.NMDS, ANOSIM, db-RDA, and heatmap plots were constructed using the vegan and pheatmap packages in R (version 4.0.2).Fungal cooccurrence networks were structured based on Spearman's correlation matrix from the OTU tables.We visualized fungal co-occurrence networks using Gephi (v 0.9.0;Mendes et al., 2014) and analyzed the topological parameters, including node number, edge number, average degree clustering coefficient, and network density.significantly lower than that in PQMF of the middle successional stage, whereas the TK content was significantly lower in QA of the later successional stage (p < .05).Comparing soil properties among non-forest, mixed forests, and R. pseudoacacia forests, the pH was the highest in the non-forest land and significantly higher than that in R. pseudoacacia forest soil.The contents of SOM and NO 3 − in the mixed forest soil were significantly higher than those in the nonforest soil but were not significantly different from those in the R.

| Soil properties
pseudoacacia soil.The TK content in the non-forest soil was significantly higher than that in the mixed forests and R. pseudoacacia forest soils.The NH 4 + content in the R. pseudoacacia forest soil was significantly lower than that in the non-forest and mixed forest soils, but the AP content was significantly higher.
The α-diversity index of fungi did not change significantly across different succession stages, but significantly changed between mixed forests and R. pseudoacacia forests (Figure 3).The Shannon, ACE, and Chao indices of non-forest, mixed forests, and R. pseudoacacia forests were significantly higher than those in three successional stage forests.The ACE and Chao indices in the mixed forests were significantly higher than those in the R. pseudoacacia forest.
The fungal β-diversity varied significantly among the three succession stages (p < .022)and between natural and planted forests (p = .039)(Figure 4 and Table 1).We further analyzed the major groups that produced differences in the fungal community structure of different forest types.Among the three succession stages, only the relative abundances of Basidiomycota and Ascomycota showed significant differences at the phylum level, whereas those of Mortierellomycota, Rozellomycota, and Mucoromycota showed no evident changes (Figure 5a).The relative abundance of Basidiomycota was significantly higher in the later successional stage QA (80.85%) soil than in the early successional stage PD (62.19%) and the middle successional stage PQMF (61.15%) soils (Table S1).Meanwhile, Ascomycetes showed the opposite pattern and were significantly lower in the QA (15.31%) than in the PD (31.82%) and PQMF (33.87%) soils.At the genus level, the relative abundance of Russula increased significantly in the QA soil during the late succession stage (Figure 5b).Between the natural mixed forest and planted R. pseudoacacia forest, Basidiomycota was significantly higher in MF (44.53%) than in RP (24.39%), whereas Ascomycota (37.11% vs. 55.57% in MF and RP, respectively) and Mucoromycota (0.40% vs. 1.20% in MF and RP, respectively) decreased significantly (Table S2).At the genus level, only Penicillium was significantly more abundant in the R. pseudoacacia forest than in the mixed forest (Figure 5b).

| Co-occurrence network complexity of fungi
Co-occurring networks of soil fungi were constructed at the OTU level to explore the effects of forest vegetation on fungal interactions (Figure 6 and Table 2).We found that the co-occurring networks of soil fungi differed with changes in the succession stage.These topological indices of the network, that is, average degree, density, and average clustering coefficient, showed a difference in the soil fungal community among the three forest succession types (Table 2).The number of nodes in the network first increased and then decreased with forest succession, reaching its highest value of 92 in PQMF in the middle stage of succession.Meanwhile, the number of edges, average degree, density, and average clustering coefficient first decreased in the middle stage of succession (PQMF) and then increased in the later stage of succession (QA).This indicates that the network complexity was lower in the middle stage of succession.
Between mixed forests and R. pseudoacacia forests, the nodes and edges of RP were higher than those of MF, but the average degree, density, and average clustering coefficient were lower.
However, the average degree, density, and average clustering coefficient of both were lower than those of non-forest.

| Functional groups of soil fungal community
The trophic modes of soil fungi in the different forest types were obtained using FUNGuild (Figure 7).The relative abundances of The relative abundance of ECM in the later successional stage QA was significantly higher than that in the early successional stage PD and the middle successional stage PQMF.In contrast, the number of saprotrophs significantly decreased in the early and middle successional stages.The relative abundances of animal pathogens, fungal parasites, endophytes, and plant pathogens were significantly higher in mixed forests and R. pseudoacacia forests than in nonforests.However, soil saprotrophs and AMF were significantly lower in mixed forests and R. pseudoacacia forests than in non-forests.
Between mixed forests and R. pseudoacacia forests, ectomycorrhiza in MF was significantly higher than that in RP, but animal pathogens and undefined saprotrophs were significantly lower.

| Correlations between fungal community and soil properties
Db-RDA showed that the fungal community structure of all 30 forest soil samples varied significantly with pH, TK, AP, and NH 4 + (Figure 8).TK content may be the major factor explaining the distinction in fungal communities among the three successional stages, whereas pH, AP, and NH 4 + may be the major factors explaining community variations between natural and planted forests.
Spearman's correlation analysis was performed separately in the three successional stages and in mixed forest and R. pseudoacacia forests to examine the relationship between soil properties and α-diversity, relative abundance of fungal taxa, and function (Figure 9).During the successional stages, only the ACE index was negatively correlated with pH (p = .045;Figure 9a).

TA B L E 2
The topological features of the fungal co-occurrence networks.
Planted forests may further affect fungal diversity by changing soil physicochemical properties such as pH, TN, NH 4 + , AP, and AK.The correlation analysis results also indicated that soil properties were significantly correlated with α-diversity (Figure 9b).
The soil fungal β-diversity showed significant changes both during the forest succession stages and between natural mixed forest and planted forest of R. pseudoacacia (RP) (p < .05).This is consistent with the finding of Jiang et al. (2021) for vegetation succession in boreal forests.Different plant species select for specific soil microbial communities (Zhang, Liu, et al., 2023).This is likely because plants secrete different root exudates, that is, amino acids, organic carbon, polysaccharides, and leteolin, during their growth, which can promote or restrict microbial communities (Lombardi et al., 2018).structure and bacterivory.However, these biotic factors were not involved in this study, and it is necessary to further consider the effects of biotic interactions on fungal communities in the future.

| Specific changes in the fungal community
We further analyzed the changes in fungi at different taxonomic levels.At the phylum level, Basidiomycota and Ascomycota were the most dominant phyla in all forest soils.This is consistent with other forest soils that typically contain two main and ubiquitous fungal groups (Dong et al., 2021;Štursová et al., 2020).They can degrade a fraction of plant residues and promote soil carbon accumulation (Clemmensen et al., 2013) contributing to an increase in Basidiomycota (Jiang et al., 2021;Yan et al., 2020).In Basidiomycota, the differences were mainly in the Russula genus, which was significantly higher in the late succession stage than in the early and middle succession stages.Russula were identified as key ECM fungi that can establish a large mycelial network to participate in nutrient cycling and energy metabolism in host plants (Twieg et al., 2007).Twieg et al. (2007) found that Russula had a high relative abundance in Q. acutissima forests, a typical tree species associated with ECM fungi.This may explain the higher Russula in QA of the late succession period in our study.
Ascomycota are able to decompose cellulose and hemicellulose in the litter (Weber et al., 2011).In this study, the relative abundance of Ascomycota in the later successional stage (QA) was significantly lower than that in the early successional stage (PD) and the middle successional stage (PQMF).Similar results were found in studies on boreal forests in China, where Ascomycota abundance was significantly lower in QA than in mixed forests (Q.acutissima and R. pseudoacacia) in November-sampled soils (Dong et al., 2021).However,  (Coutinho et al., 2012;Efthymiou et al., 2018).The higher relative abundance of Penicillium in the planted forest RF may explain the significant increase in soil AP content (Figure 1).Owing to the low bioavailability of phosphorus, a shortage of AP in soils often becomes a key factor restricting plant growth and development (Vitousek et al., 2010).Therefore, the enrichment of Penicillium in R. pseudoacacia forests is beneficial for promoting the availability of soil nutrients, which may play an important role in maintaining R. pseudoacacia forest growth and health.

| Functional changes in the soil fungal community
We found that ECM increased in the late succession stages (QA).
The main reason for this phenomenon may be that the QA forest in the late succession stage is a typical tree species associated with ECM fungi, forming a large mycelial network with ECM fungi that participate in nutrient cycling (Tedersoo et al., 2010).In contrast, soil nutrient contents decreased significantly in the late succession period, such as TK, NO 3 − , and AP.This may provide an oligotrophic environment favorable for the colonization of ECM fungi (Khalid et al., 2021).ECM fungi can promote water and nutrient absorption in host plants through the mycelial network, enhance the stress tolerance of host plants, and facilitate community succession (Kałucka & Jagodziński, 2016).In contrast to ECM fungi, the relative abundance of saprotrophs decreased with succession.
Saprotrophic fungi can accelerate the decomposition of SOM to facilitate nutrient cycling in terrestrial ecosystems, among which plant saprotrophs are mainly involved in plant decomposition (Bello et al., 2021).Increasing soil carbon sources can stimulate saprotroph fungal growth, which may result in an increase in the relative abundance of saprotrophs (Marañón-Jiménez et al., 2021).
In this study, the increase in the relative abundance of plant saprotrophs in the early successional stage may be attributed to the higher carbon content of the leaf litter in coniferous P. densiflora forest (PD) than in broad-leaved Q. acutissima forest (QA) (Lu et al., 2022).This indicates that plant species play an important role in regulating the function of soil fungal communities during forest succession.
The relative abundance of animal pathogens and undefined saprotrophs was higher in R. pseudoacacia forest soil than in mixed forest soil, whereas ectomycorrhizal abundance showed the opposite trend.This significant change in fungal function may be attributed to the elimination of major habitats for plant symbionts, such as ECM, following the removal of the forest floor (Tedersoo et al., 2003).
ECM and saprobic fungi are most sensitive to logging disturbances (Hartmann et al., 2012).In addition, R. pseudoacacia is symbiotic with AMF (Yang et al., 2015).Arbuscular mycorrhizal trees have a higher relative abundance of saprotrophs and pathogens (Eagar et al., 2022;Fang et al., 2020).More attention should be paid to changes in pathogenic fungi under planted forest management.

| CON CLUS ION
Our study demonstrates that the development of different forest types via natural and planted restoration processes exerts a profound impact on soil fungal diversity, community composition, and functional group dynamics.Forest succession was observed to affect the fungal community composition, yet it did not alter the α-diversity.In this typical succession sequence from coniferous P. densiflora to broad-leaved Q. acutissima forests, a notable increase in ECM and a decrease in saprophytic fungi were evident during the later successional stages.Conversely, reforestations were found to reduce overall soil fungal diversity and alter community composition.In terms of functional groups, R. pseudoacacia forests decreased the relative abundance of ECM and increased the number of animal pathogens.It is worth exploring the role of specific differences in fungal taxa in forest succession and reforestations in future research.
gal communities to reforestation to promote the healthy development of forest communities and maintain forest ecosystem stability.Kunyu Mountain, located east of the Jiaodong Peninsula, China, is the native land of Pinus densiflora Siebold & Zucc.and is the largest and best-protected natural distribution center of P. densiflora worldwide.This area plays a vital role in water and soil resource conservation and is the most important ecological screening on China's Shandong Peninsula.The forest community mainly comprises P. densiflora, interspersed with mixed forests of P. densiflora and Quercus acutissima Carruth., as well as stands of Q. acutissima.This distribution exemplifies a typical successional sequence transitioning from coniferous to broad-leaved forests.The area also features natural mixed forests, which include species such as pine, oak, and Robinia pseudoacacia L., alongside planted pure forests of R. pseudoacacia.The changes in plant species under forest secondary succession and reforestation may have significant effects on the T A X O N O M Y C L A S S I F I C A T I O N Biodiversity ecology, Community ecology, Microbial ecology, Soil ecology | 3 of 14 ZHU et al.
Forest vegetation succession significantly affected the contents of TN, TK, NO 3 − , and AP, whereas pH, SOM, TP, NH 4 + , and AK showed no significant changes.NO 3 − and AP contents increased significantly and then decreased during forest succession (Figure 1), reaching the highest levels in the PQMF during the middle period of succession.The TN content in PD of the early successional stage was F I G U R E 1 Comparison of soil properties among different forest types on Kunyu Mountain.Lowercase letters indicate statistically significant differences.MF, deciduous broadleaved mixed forest; NF, non-forested land; PD, P. densiflora; PQMF, pine-broadleaf mixed forest of P. densiflora and Q. acutissima; QA, Q. acutissima; RP, R. pseudoacacia.

F
I G U R E 5 Comparisons of the relative abundances of the major taxa of fungi in different forest types at the phylum (a) and genus (b) levels.Lowercase letters indicate statistically significant differences.MF, deciduous broadleaved mixed forest; NF, non-forested land; PD, P. densiflora; PQMF, pine-broadleaf mixed forest of P. densiflora and Q. acutissima; QA, Q. acutissima; RP, R. pseudoacacia.
010) and Penicillium (p = .010)were negatively correlated with NH 4 + , while Mortierella was negatively correlated with pH (p = .048).Functionally, ectomycorrhizal was positively correlated with NH 4 + (p = .027)and negatively correlated with AP (p = .019),while animal pathogens were positively correlated with AP (p = .002).The plant pathogen was negatively correlated with pH (p = .008),but positively correlated with SOM (p = .002)and AP (p = .040). 4 | DISCUSS ION 4.1 | Different responses of α-and β -diversities to forest type Forest type significantly affected soil fungal diversity in this study, but α-and β-diversity showed different patterns.Forest succession did not significantly affect the α-diversity of fungi.Soil fungal α-diversity has been shown to be significantly affected by edaphic variables in forest ecosystems, especially soil pH and nutrient content (Vasco-Palacios et al., 2020; Wang et al., 2015).In this study, nonsignificant changes in fungal α-diversity in the three succession stages may be related to there being relatively little change in soil properties, namely, pH, SOM, TP, NH 4 + , and AK (p > .05).The results of the correlation analysis also showed no significant influence of soil physicochemical properties on fungal α-diversity during forest succession (Figure 9a).However, between natural mixed forest and F I G U R E 7 Changes in soil fungal functional groups among different stages of succession and between natural and planted forests.Lowercase letters indicate statistically significant differences.F I G U R E 8 db-RDA plots for fungal assemblages showing the covering relationship between soil properties and community structure.MF, deciduous broadleaved mixed forest; NF, nonforested land; PD, P. densiflora; PQMF, pine-broadleaf mixed forest of P. densiflora and Q. acutissima; QA, Q. acutissima; RP, R. pseudoacacia.
Li et al. (2022) found that Populus secretes strigol and phenolic compounds, which are key drivers of bacterial and fungal communities in the rhizosphere.These findings indicate that changes in the dominant plant species have a strong influence on the fungal community.Moreover, the changes of environmental factors also significantly affected the community structure of fungi.The RDA plot showed that TK, pH, AP, and NH 4 + were important driving factors affecting the community structure of soil fungi in different forest types.However, it is worth noting that only a small fraction (about 18%) of variation in fungal community structure was explained by the selected environmental variable.In fact, besides the environmental factors, biotic interactions are known as a major control on microbial diversity (Zhou & Ning, 2017), such as predation, competition, allelopathy, parasitism, and mutualism.Li et al. (2023) found that the bacterial diversity indices and community structure were much more accurately explained by a number of eukaryotic characteristics than by the measured environmental variables and/or spatial variables alone.Over 70% of total variation in bacterial community structure could be explained by three categories of biotic interactions: parasitism, fungi-bacterial competition, and trophic F I G U R E 9 Spearman's correlations between soil properties and fungal α diversity, dominant genera, and functional groups at three successional stages (a) and in natural and planted forests (b).Asterisks and double asterisks indicate p values <.05 and <.01, respectively.The values of the correlation coefficients are indicated according to the color bar.
(Torres-Cruz et al., 2018)fferences occurred in the Penicillium genus, which is a highly competitive species because it can outcompete other microbes by producing large numbers of asexual spores and surviving in disturbed environments(Torres-Cruz et al., 2018).This may account for the relatively high abundance of Penicillium in the planted forest soil used in the present study.Penicillium is typically a phosphate-solubilizing fungus, and many strains can solubilize forms of phosphate that are unavailable to plants (Chen et al., 2017)d forest and R. pseudoacacia forest, we found that Basidiomycota decreased in RP, while Ascomycota and Mucoromycota increased.Ascomycota are generally abundant in environments with high disturbance or strong human intervention compared to virgin forest soils(Chen et al., 2017).